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The Affleck-Dine mechanism, which is one of the most attractive candidates for the baryogenesis 
in supersymmetric theories, often predicts the existence of baryonic Q balls in the early universe. 
In this scenario, there is a possibility to explain the observed baryon-to-dark matter ratio because 
Q balls decay into supersymmetric particles as well as into quarks. If the gravitino mass is small 
compared to the typical interaction energy, the longitudinal component of the gravitino behaves like 
the massless goldstino. We numerically calculate the goldstino production rates from Q balls in the 
leading semi-classical approximation without using large radius limit or effective coupling. We also 
calculate the quark production rates from Q balls in the Yukawa theory with a massive fermion. In 
deriving the decay rate we also take into account the scalar field configuration of the Q ball. These 
results are applied to a realistic model in the gauge-mediated supersymmetry breaking and yield the 
Oh branching ratio of the Q ball decay into the gravitino. We obtain the branching ratio much smaller 

, than the one estimated in the previous analysis. 

m 

CS| . I. INTRODUCTION 

In cosmology, the origin of dark matter and the baryon asymmetry is very important problem which can not be 
[ explained by the Standard Model of elementary particle physics. Supersymmetric (SUSY) extensions of the Standard 
JL . Model can explain this problem. In minimal supersymmetric standard model (MSSM) the lightest SUSY particle 
" (LSP) is stable and is natural candidate for the dark matter. Besides, the baryon asymmetry can be produced by the 
Affleck-Dine mechanism [l], The Affleck-Dine mechanism produces a scalar field condensate with baryon number. 
In many models such as the gauge-mediated SUSY breaking (GMSB) models, this condensate feels spatial instabilities 
| | and fragments into nontopological solitons, Q balls 0-0] ■ 

In the GMSB, the LSP is gravitino. If the charge of the Q ball is large enough, the Q ball is stable. Since the 
y—{ ' stable Q ball has astrophysical problems [1, 0, we do not consider this case. On the other hand, if the charge of the 
Q ball is small enough, the Q ball decay into hadrons and gravitinos. If Q balls decay only via squark — > quark + 
r**"* ' gravitino [ToL [ll[ , the produced gravitino number density is equal to the produced quark number density by R-parity 
*n , conservation. Then the gravitino mass should be ~ 1.7GeV to explain the observed baryon-to-dark matter ratio 
~ 1/5. However, it is pointed out that the process as squark + squark — > quark + quark occurs via heavy gluino 
exchange and is main decay mode if decay into a squark is kinematically forbidden [12| . In this case, it becomes 
important to calculate the branching ratios into the gravitino and quarks to estimate the baryon-to-dark matter ratio. 
In this article, therefore, we derive the production rates of gravitinos and quarks from the Q ball in the GMSB 
. model (the most interesting case) and in the gravity-mediated SUSY breaking model. The Q ball decay was first 
studied by Cohen et al. [13], who considered the Yukawa theory and treated the scalar field as the classical Q 
ball background field in the leading semi-classical approximation and calculated the fermion production rate by the 
Bogoliubov transformation between creation and annihilation operators at t ±oo [13| . In the case of gravitino, if 
the gravitino mass is small compared to the typical interaction energy, the longitudinal component of the gravitino 
behaves like the massless goldstino and has the derivative interactions with chiral multiplets. Therefore we consider 
fermion fields (one of which is a goldstino) with derivative coupling to a scalar field. Once we treat the scalar field 
as the classical background field, the interaction terms become at most second order for the fermion fields and this 
system can be solved numerically in the same way as Cohen et al. did. The gravitino production rate from Q balls 
was first calculated in Ref. |12], where the effective coupling estimated from the decay rate of squark — > quark + 
gravitino and large radius limit are used. We derive the gravitino production rate without relying on the effective 
coupling and large radius limit. We also consider two fermion fields one of which has a mass term in the Yukawa 
theory. This system is almost equivalent to the one considered in Refs. [g, The interaction of Q balls with 

ordinary matter was studied in Ref. Q , which showed that quarks are reflected as antiquarks with a probability ~ 1 
but did not refer to the Q ball decay. In Ref. [l4| , the authors estimated the quark production rate using the effective 
interaction after integrating out the heavy particle. However, in this paper we calculate the rate of the Q ball decay 
into quarks without integrating out the heavy particle. Furthermore, in deriving the decay rate we also take the scalar 
field configuration of the Q ball into account. In most of the previous studies, a step function is used as the scalar 
field configuration. We carefully examine the differences between the realistic and step-function configurations. 
This paper is organized as follows. In SecHH we briefly review the property of Q balls in the GMSB and the gravity 



X 



2 



mediated SUSY breaking models. In Sec. IIII1 we review the method to calculate the decay rate of Q balls and discuss 
the differences among the decay rates for some Q ball configurations in the Yukawa theory. We apply the method to 
calculate the goldstino production rate from the Q balls in the GMSB and in the gravity mediated SUSY breaking 
models in Sec. IIVI In Sec. [Vj we also calculate the Q ball decay rate in a theory with a massive fermion. Then, we 
apply our results to gravitino and quark production from Q balls in SecEU Sec. IVIIl is devoted to the conclusion. 

II. Q BALL SOLUTIONS 

In this section, we consider a complex scalar field theory with a global U(l) symmetry. The Lagrangian density is 
written as 

C=-d ll ct>*d^(j>-V{(t>). (1) 
The conserved Z7(l) charge and the energy are given by 

Q = iJ(<f>*do<P-(f>do<P*)d 3 x, (2) 

E = J (do^do^ + d^d^ + Vi^^x, (3) 

respectively. The scalar field configuration which minimizes the energy at a fixed charge Q is obtained by minimizing 

E + uJQ-i [ (^d cf> - cj>d Q ct>*) d 3 x) , (4) 



where ujq is a Lagrange multiplier. We can immediately determine the time dependence of the scalar field configuration 
and obtain 4>(r,t) — 0(r)e _Mo *. Then, taking a spherically symmetric ansatz <p(r) = <p(r), we get the radial part of 
the configuration by solving the following equation: 

with the boundary condition 0'(O) = 0. Here, if ui^ < min^ (r|p) < V"(0)/2, we get the spatially localized 

configuration, called Q ball (l5j . 

In MSSM, there are many flat directions in the scalar potential. The flat directions are combinations of squarks, 
sleptons and Higgs, but they are lifted by the SUSY breaking effect. In the early universe, the Affleck-Dine mechanism 
produces a scalar field condensate with baryon number using one of these flat directions. If a Q ball solution exists 
for this flat direction, the scalar field condensate fragments into Q balls 16]. In the following subsections, we review 
the properties of the Q ball in the GMSB and the gravity-mediated supersymmetry breaking models. 



A. Properties of the Q ball in gauge mediation 



In the GMSB model, the gauge fields acquire large masses for g<j> ^> M m , where M m is the messenger scale, and g 
generically stands for the standard model gauge coupling. Then the transmission of SUSY breaking effect is suppressed 
and the flat directions flatten out for gcj) ^> M m as [l7| 



V = 



2mlMl 
9 2 



log 1 



(6) 



where m s is the soft mass scale. In the GMSB gravitino mass ro 3 / 2 is much smaller than m s . In this potential, there 
always exists a Q ball solution. The mass of Q ball Mq, the typical size of the Q ball R, ujq, and the field value at 
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the center of the Q ball 4>q = <f>(r — 0) are written as 



AJ2ttc fmMn x 1/2 



M Q c ^ i^Hj Q3/4 ; (?) 
^M(^T-Y /2 Q^U-), (8) 



w 



where the parameter c is fitted as \U 



2c \m s M m J \ ujq 



2=7- ( J V2 £T 1/4 , (9) 

c/^y /2 gl/4 / c!! a\ 

2tt V 5 / V 5 w o / 



4.81og(m s /a;o) + 7.4. (11) 



The scalar filed configuration is well approximated by <j> = 4>q sm ^° r ^ for r < R. However, this configuration is not 
smooth at r = R [see Eq. ©]. This cause some difficulty in deriving decay rate of the Q ball into goldstinos. Thus 
in order to obtain the smooth configuration we solve Eqs. ([5]) and © numerically. We show the numerical solution 
in Fig. [TJ The second derivative drastically changes around rujQ ~ tt, but it is finite and nonsingular, as expected. 




FIG. 1: Numerical solutions of Eq. pa). cf>(r), <j>'{r) and (j>"(r), in the GMSB with m a /ojo = 10 3 . The vertical axis is normalized 
by (cm a M m / guJo) for <j>(r) and c/>'(r), and (cm s M m /guio x 1/10) for cj>"(r). The horizontal axis denotes the radius normalized 
by l/uo- 



B. Properties of the Q ball in gravity mediation 



In the gravity-mediated SUSY breaking model, the fiat directions are lifted by the SUSY breaking effect as 

V = ml\^(l + K\og^.y (12) 

where m s is the soft mass scale and Mq is the reduced Planck mass. In the gravity-mediated SUSY breaking gravitino 
mass is the same order of m s . The second term in the parenthesis comes from radiative correction, and typically 
\K\ ~0. 01-0.1. It is known that K < for the flat directions of the first and second generation squarks, and then 

there exists a Q ball solution The scalar field configuration is well approximated by </>(r) ~ (f)o exp ( — J, and 
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Mq, R, wo and </>o are written as 

M Q ~ m s g, (13) 

(14) 

uj ~m a , (15) 

O ^(27T 3 / 2 )- 1 /2|X|3/4 ms Ql/2 ; (16) 

III. Q BALL DECAY IN THE YUKAWA THEORY 

In this section, we consider the Yukawa theory in the Q ball background field. In this system, Cohen et al. proved 
that the Q ball decay into fermions in the leading semi-classical approximation. Here, we first review the method 
presented by them and then apply to some cases without using large radius limit. 

The Lagrangian density is written as 

C = rfifdux + riUa^v ~ {94>*XV + h.c.) , (17) 

where x an d V are two-component Weyl spinors, and cr M = (1,—a 1 ), where a l, s are the Pauli matrices. We assign 
the global U(l) charge for 0, \ an d V as 1, 1 and 0, respectively, and treat the scalar field as the Q ball classical 
background field <j> — <f>{r)e^ tulot in the leading semi-classical approximation. The discussion below is also correct for 
more general interactions like the one used in the following sections. 

When we treat the scalar field as the classical background field, this system is invariant under simultaneous time 
translations and U(l) rotations as 

l t " t + 5t \ t (18) 

The associated conserved current is 

f = x Ua» (do + ^wo) X + r]Ua^doV - 5$C, 

where is the energy momentum tensor for the fermions and is the U(l) current. Then, the conserved charge is 
not the energy Ef itself but Ef — u>qQ{, where Qf is the U (1) charge for the fermions. Using this conservation law, we 
can show that the fermion fields are created from the Q ball surface in the leading semi-classical approximation. We 
have treated the scalar field as the Q ball classical background field and the Lagrangian density is at most second order 
for the fermion fields, so we can easily solve this quantized fermionic system. Since the conserved charge is not Ef but 
Ef — u>oQf, positive and negative modes may mix with each other in this system. Thus, even if an incoming wave is a 
positive mode (annihilation operator) at t = — oo, the outgoing wave can be a negative mode (creation operator) at 
t = oo. We can regard this as the Bogoliubov transformation between annihilation and creation operators at t — > ±oo, 
and then we can calculate the non-zero number density of outgoing waves. 

We need the free field mode expansions of the fermions to consider the scattering problem of the fermion fields in 
the Q ball background field. Far away from the Q ball surface, the background field is <p — 0j and then the equation 
of motion for the mode of x oc e~ lk+t is 

(k+ - ia ■ V) x = 0. (20) 

In this article, we treat Xa, X^°S Va and rf a as column vectors and neglect the indices. Due to the rotational invariance 
of this system, we can expand the solutions by the following Pauli spinors: 



$(j,m,/=j + l/2) = 



$(j,m,/'=. 7 -l/2) 



V / j+l_-m ym-l/2 
\/2(7+T) ' 
Vj+l+m ym-H/2 

Jj+m v m — l/2\ 

Jj — m y-m+1/2 
V2j I' 



(21) 
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where Y t m are the spherical harmonics. Defining uW (i=l, 2) as 

«W(fc,i,m;r) = A^^fcrJ^m.j'J + ftW^^m,/)], (22) 

where ftj 1 ^ and ftj 2 ^ are the spherical Hankel functions of the first and second kinds, respectively, we have 
(k + ia ■ V) ttW(fe) = 0. So we obtain the solutions outside the Q ball as 



poo 

X = Y, dk+{a in (k + ,j,m)e- ik + t u ( - 1 \-k+,j,m;r) 

A „ JO 



a ut(fc+, J, m)e lk+f u' 2 '(-i; + , j, m; r) + (terms of antiparticles)}, (23) 



where a; n and a out are expansion coefficients. We expand r\ in the same way as x but denote expansion coefficients 
as c instead of a. For later use, we write the expansion of = rf a as 

POO 

J = i°2 (Va)* = E / *+{cL(fc+,i.-m)e ife +*(-l) m - U ( 2 )(-fc + ,j,m;r) 

+ cl ut (k+,j, -TO)e ife +*(-f) m -u (1) (-fc + , j, m;r) + (terms of antiparticles)}, (24) 

where we have used ia2U^' 2 \k, j, m; r)* = (— f )" l +u' 2 ' 1 ^(fc, j, — to; r), where m± = to ± 1/2. 

When we quantize the fields, a; n and a ou t become the annihilation operators for incoming and outgoing waves, 
respectively. The operators obey the Heisenberg equations of motion and are scattered by the Q ball at the origin. 
We shall consider this scattering problem. 

First, let us consider the case without Q balls. Since the solution has to be nonsingular at the origin, the radial part 
of the solution is written by the spherical Bessel functions, ji(kr), as x « u(—k,j,m;r), where u(k,j,m;r) is defined 
by substituting the spherical Bessel functions for the spherical Hankel functions in the definition of «W. Using this, 
along with ji = (h^ + h^)/2, we obtain a out = a ln . This means that the incoming wave reflects off the origin and 
becomes the outgoing wave. Of course, the number density of the outgoing wave is (0i n |aQ Ut a out |0i n ) = when the 
vacuum, |0; n ), is defined as the state of no incoming wave, i.e. ai n |0i n ) = 0. 

Next we consider the case with the Q ball background field whose center is located at the origin of the coordinate. 
The Heisenberg equations of motion are 

(25) 

jxfdrf - g(f>* X = 0. 

Since this is the combination of the linear differential equations, we can solve these equations with appropriate 
boundary condition, and then outgoing waves, a ou t and c ou t, can be written by incoming waves, a; n and c; n , like the 
previous example. However, because the background field depends on the time as <f) °c e~ lUot , the conserved quantity 
is not Ef but Ef — WoQf- Thus, the modes which mix with each other are 

f x K e -i«t = e -iM, 

]y cx e -*("-"o)t = e *fc-t. ^ > 

We are interested in situations where positive and negative modes mix with each other, because we derive the fermion 
production rate through the Bogoliubov transformation between creation and annihilation operators. Thus, we restrict 
our attention to < to < uj - 

Outside the Q ball, the fields satisfy the free equations of motion. The angular momentum conservation implies 
that the terms which mix with each other can be written by 

ia in (k + J, m )e^ k + t u^ 1 \-k + J,m;r) +a out (k + J,m)e- lk + t u^(-k + ,j,m;r), 

{c^fc_,j;-TO)e* fc -'(-ir- M ( 2 )(-fc_,j,TO,^^^ 1 } 

where a out and CQ Ut can be written by a; n and c\ n by matching the interior and exterior solutions. Using the super- 
position principle, we can write 



a out (k + ,j,m) \ fR x (k + ,j) T x (k+,j)\ ( a in (k + ,j,m) 
(-ir-ct ut (k-,j,-m)J \T„(k-J) Rr,(k-,j)J {(-ir-cl(k-,j,-m))> 



(28) 
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where Ri and Ti are coefficients fixed later by matching the interior and exterior solutions and do not depend on 
m due to rotational invariance. We can regard this as the Bogoliubov transformation between t — > ±00. The 
anticommutation relations {ai n ,aj n } = {a ou t,aJ ut } = {ci n , cj n } = {c out ,cl ut } imply that this translation matrix is a 
unitary matrix. Especially, we have 

(\T x (k + ,j)\ 2 = \T v (k-,j)\ 2 , 

\\R i (k + ,j)\ 2 + \T i (k.,j)\ 2 = l, (i = X ,v)- [ ' 

We define the vacuum, |0i n ), as ai n |0i n ) = Cj n |0j n ) = at r — > 00, where the incoming waves are not affected by the 
Q ball. The incoming waves are scattered by the Q ball and this scattering process is described by the Heisenberg 
equation of motion. Then, we get Eq. (|28p and the number density of the outgoing x waves is 



(O in |aL t (fc+,J, m)a mt (k' + ,j',m')\O in ) = (-l) m - +m '-r;(fc + , jOT x (fc;,/)(O in |c in (fc_, j, -m)^^,/, -m')|0 in ), 

= l^x( fc +)l 2 ^( fc + ~ k' + )6jj'5 mim >. 

This proves that the outgoing x waves are created by the presence of the Q ball background field. Similarly, the 
outgoing rj waves are also created. The above unitarity condition, |T x (fc + , j)| 2 = |T^(fc_,j)| 2 , of Eq. (|29|) implies that 
the process can be represented by <I>bg(E — w o) — > x(E = k+) + v{E = k- = uj — fc+), where <Pbg is the Q ball 
background field, and E is the energy of each field. The second condition, \R x (k+, j)\ 2 + \T x (k-, j)\ 2 — 1, of Eq. (|29|) 
implies that the production rate is bounded above due to the Pauli exclusion principle. Summing over the states and 
using (5(0) = T/2-7T, we obtain the following production rates of the fields from the Q ball: 

J t N *=J2 2-(2i + l)m(fc,j)| 2 , (i = X ,v)- (31) 

■ .1 Q 



3=1/2 



Finally, we need to determine the coefficients Tj by matching the interior and exterior solutions. Inside the Q ball, 
the Heisenberg equations of motion are given by Eq. (j44|) . With use of Eq. (f26|) . the equations are rewritten as 



' (k+ - ia ■ V) x ~ 9<l>(rW = °> 
(-A;_ + ia ■ V) 77 f - g<j){r)x = 0, 
The conservation of the angular momentum implies that the solutions can be expanded as [if 

fx = fx ( r ) $ U: m ' l ') + i9x ( r ) $ U: m > » 

I = fv ( r ) * (i. TO > + Wv ( r ) $ C?> TO ' ■ 



(32) 



(33) 



Then, using 



a.V<I»l./.»/..,'±J/2) = #(./.m../ ! 1/2) (|:+ 1±0 r +1/2) )- <-H) 



we obtain the four first order differential equations. Since the solutions have to be nonsingular at r = 0, the boundary 
conditions are f[ (r = 0) = g\ (r = 0) = for i — x and rj. We can get two independent solutions numerically for the 
given scalar field configuration (f>(r). Then, matching the interior and exterior solutions and using Eqs. (f2"Tl) and (|28l) . 
we obtain the coefficients Ti. The matching condition is simply setting the solutions equal to each other at sufficiently 
large r where 4>(r) ~ is satisfied. 

A. Q ball decay rate for Rojq — s> 00 in the Yukawa theory 

Using the above technique and the large radius limit (Rluq — > 00), Cohen et al. derived the Q ball decay rate for 
the Yukawa interaction [l3| . The scalar field configuration was taken to be the step function as 

<t>tr)=mR-r)=\^ °< r * R > (35) 

0, R < r. 
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In the limit of g<f>a/u)Q 3> 1, the production rates are saturated by the Pauli exclusion principle and written as 

d Ni= (dN\ s ^ for ^ >>1; (36) 



dt \ dt / sat 247r UjQ 
where i = \i V- On the other hand, in the limit of g<fio/ujQ "C 1, they showed that the production rates are 

* Ni *to(°*)x(*pl , far < 1, A* oo. (37) 
dt \ uj J \dt ) sat w 

We can see the physical meaning of this behavior in the following way. The penetration length of incoming waves 
inside the Q ball is ~ l/(g0o) [l3| , so the effective volume of the interaction near the surface of the Q ball is 
V e s ~ 47ri? 2 /(50o). Then, the Q ball decay rate, which is roughly the decay rate of the scalar field (~ g 2 (^o) times 
the charge density times the effective volume, is ~ {g 2 ^o) x (^a^o) x ^cff • 



B. Q ball decay rate for Rljq ~ 1 in the Yukawa theory 



In the realistic models of the previous section, the radius of the Q ball R is ~ 1/wo, and the effective volume of 
interaction become the whole region of the Q ball, V e s ~ AttR 3 /3 14]. The parameter dependences of the Q ball 
decay rate are, therefore, quite different from the previous subsection. We derive a fitting formula for the decay rate 
in the limit of g<f>o/ujQ *C 1 and Rluq ~ 1. The decay rate are shown in Figs. [5]and[3]as a function of R and gifio/coo, 
respectively, and can be fitted as 

( d -§) cl2f^) 2 (i?w -1.9)xff ) , for^«l, R. >2, (38) 

V dt J stcp function \oj J \dt J sat ujq 

for the step-function type of scalar field configuration. 
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FIG. 2: Production rate from the step-function type of Q ball as a function of Ruo with g<j>o/ijJo = 0.01 in the Yukawa theory. 
The vertical axis is normalized by the saturated rate of Eq. (|36p . This result can be fitted as (dNi/dt) ~ 1.2 x 10~ 3 (7?cj - 
1.9) x (dN / dt) sa x for Ruiq > 2. The production rates of i — \ and rj are the same. 

We numerically calculate the fermion production rates for not only the step-function type of Q ball configuration 
but also the two types of Q ball configuration introduced in Sec. HU the numerical configuration shown in Fig. [T] for 
the GMSB model and (f>(r) — <$>oe~ r l' 2R for the gravity mediated SUSY breaking model. The results are shown in 
Fig. [3] as a function of gfo/ojo with Rujq = n. We can explain the differences among the production rates for the 
each type of scalar field configuration. In the limit of g(f>o/^o < 1, Q ball charges should be taken into consideration 
because Q ball decay can be regarded as a collection of decay of the scalar field. Since the charge density is given by 
2cjo0o a t the center of the Q ball, the Q ball charge is roughly estimated as (87r/3)i? 3 wo0o- However, for the realistic 
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configuration of the Q ball, this estimation is not a good approximation. The ratio of the actual total charge Q to 
the rough estimation is 



3Q 



8nuj R 3 (t)l 



fl, 



3 

3^ 



for step function, 



— -it, for gauge mediation, 

2tt z 



for gravity mediation, 



(39) 



for Rujq = 7r. Thus, we can approximate the production rates for each type of configuration as 



dt 



jL x ( dNt 

2tt 2 X V dt 

30F ^ fd_N, 
~dt 



step function 



step function 



for gauge mediation, 
for gravity mediation, 



(40) 



for gtfio/uJo *C 1 and Rloq = -k. On the other hand, for g^o/coo 3> 1, the Q ball decay rates are saturated by the Pauli 
exclusion principle. Since the penetration length of incoming waves inside the Q ball is ~ l/(g4>), the Q ball decay 
rates depend on R' , where R' is determined by g4*(R')/uJo ~ 1. Thus, from Fig. [TJ the decay rate for the gauge- 
mediation type of Q ball becomes the same as that for the step- function type of Q ball in the limit of gfio/ujQ 1. 
Since the gravity-mediation type of configuration is 4>(r) — 4>oe~ r / 2R , the effective radius logarithmically increases 
as R' ~ R(2 \og(gcj)o/uJo)) 1 ^ 2 and the decay rate for that of Q ball also increases as (dN/dt) sa .t\R-^R' — 2 log(g0 o / w o) x 
(dN/dt) sat in the limit of gcf>o/uJo ^S> 1. These considerations explain the decay rate of Fig. [3l The agreement on the 
decay rate with the step-function type of Q ball is very good for the gauge-mediation type of one. On the other hand, 
the agreement is not so good for the gravity-mediation type of Q ball. This disagreement comes from the fact that 
the step-function type and gauge-mediation type of Q ball are thin wall configurations, but the gravity-mediation 
type of one is a thick wall configuration. We conclude that the decay rate for the gauge-mediation type of Q ball can 
be approximated as 



dNi 
dt 



gauge mediation 




step function 



for < 1, 

UJQ 

for 3> 1, 

UJQ 



(41) 



where (dN/dt) sat is given by Eq. (|3T>1) . and (dNi/dt) s t cp function can be fitted as Eq. (f3"5)) . The decay rate for the 
gravity-mediation type of Q ball can be fitted as 



dNj 
dt 



3 x 



dNi 
dt 



step function 



gravity mediation 



5.9 + 1.75 log ( ) - 0.02 



log m 

UJQ 



dN 
~dt 



for < l, 
for ^ » 1, 



(42) 



for Rloq = 7r from our numerical calculation. 



IV. Q BALL DECAY RATE INTO GOLDSTINOS 



We apply the method reviewed in the previous section to derive the goldstino production rate. When global SUSY 
is spontaneously broken, there is a goldstino and its interaction with a chiral multiplet at low energy is given by 

Ant = ~ (vd, (K^X) d^*) + 9 M (d^x^V) 7?t) , (43) 

where (F) is the SUSY breaking F term, and cr M = (ljU 1 ). Here, 77, \ an d 4> are the goldstino, the chiral fermion and 
its superpartner, respectively. We assign the global U(l) charge for (f),x an d T) such as 1, 1 and 0, respectively, and 
treat the scalar field as the background field = (f){r)e~ lulat . 
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step function 
gauge mediation 
gravity mediation 



0.01 0.1 1 10 100 



FIG. 3: Production rates from Q ball as a function of g<j>o/ujo with Rluq = n in the Yukawa theory. The types of scalar field 
configurations are taken to be the step-function type (red line), the gauge-mediation type (green dashed line) and the gravity- 
mediation type (blue dotted line). The vertical axis is normalized by the saturated rate of Eq. (|36[) with Rujq = n. The ratio 
of the step-function type of Q ball to the gauge-mediation type of one to the gravity-mediation type of one on the production 
rate is 1 to 1/7 to 3 for g(f>o/ijJo <C 1. The production rate for the gravity-mediation type of Q ball can be approximated as 
(dNi/dt) ~ (5.9 + 1.75 log (g4>o/cuo) — O.O2(log(<;0o/^o)) 2 ) x (dN/dt) sat for gcj)o/uJo 3> 1. The production rates of i = x an d V 
are the same. 



The analysis is the same with the previous section once we replace the Heisenberg equations of motion of Eq. (j2"5 
with 



S^ctVgW = 0, 



2 1 



(44) 



With use of Eq. (121)1) , this equations are rewritten as 
'(k+ - ia ■ V)x 

+ — — (w #- ~ 20' d r + w <t> (i<J ■ V) + 0'fc_ (ia r ) - 4> (ia r ) (ia ■ V)) = 0, 
(-k- +ia ■ 

+ L>l<j> + i^- (rV) - CJ #+ + u c/) (zcr • V) + fc+0' (ia r ) - (j)' (icr r ) (in ■ V)\ X = 0, 



(45) 



where a r = a ■ r/r and (/>' = d<fi(r)/dr. These equations correspond to Eq. (1521 in the previous section. We need the 
following relation, in addition to Eq. (|34[) . to simplify the equations: 

(46) 



<5v$ (j, m,j ± 1/2) = $ (j, m, j T 1/2) . 



In this section, we use the two types of the Q ball configuration 4>(r) introduced in SecHU the numerical configuration 
shown in Fig. [1] for the GMSB model and <p(r) = (t>oe~ r / 2R for the gravity mediated SUSY breaking model. 



A. Goldstino production rate from the Q ball in gauge mediation 



In the GMSB model, the scalar field configuration is taken as the numerical configuration shown in Fig. [T] We show 
an example of the energy spectrum of outgoing \ waves in Fig. [4] Using |T x (fc_, j)| 2 = \T v (k +7 j)\ 2 of Eq. (|29|) . along 
with fc_ = luq — k + , we can also see the 77 production rate from this figure. The energy dependence of the production 
rate is slightly asymmetric by the replacement of k + — > loq — k + , so either x or V gets more energy from Q ball than 
the other. 
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FIG. 4: Energy spectrums of the outgoing \ waves \T x {k+, j)\ 2 for (f>wa/{F) = 1 (red line) and ip^o/iF) = 10 (green dashed 
line) with j = 1/2 in the GMSB. 




§q(DqI<F> 

FIG. 5: Production rate from the gauge-mediation type of Q ball interacting with the goldstino as a function of 4>qujo/(F). 
The vertical axis is normalized by the saturated rate of Eq. (|36|) with Ru)o = n. The production rates of i = % an< l V are the 
same. 



The production rate is shown in Fig. [5] as a function of 4>o^o / (F) and it can be written as 



d M \"\(F) J \dtj s&t ' - (F) <lf 

where (dN/dt) sat is given by Eq. (l36l) . We can also get this behavior using Eqs. (j38| and (|4Tj) if we naively estimate 
derivatives in the interaction term as d/dr ~ 1/R and <9/<9i ~ ujq and use i?u;o ~ ir. 



B. Goldstino production rate from the Q ball in gravity mediation 

In the gravity-mediated SUSY breaking model, the scalar field configuration is taken as <fi(r) = (j> exp (— r 2 /2i? 2 ), 
where R — |_K"| _1 / 2 u;(^ . We show an example of the energy spectrum of outgoing \ waves in Fig. [5] The production 



11 



rate is saturated for k + ~ wo/2 by the Pauli exclusion principle. Using \T x (k-, j)| 2 = \T v (k + , j)\ 2 of Eq. (|29)) . along 
with fc_ = ujq — k + , we can also see the 77 production rate from this figure. 
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FIG. 6: Energy spectrums of the outgoing x waves \T x (k+, j)\ 2 for 4>ujo/{F) = 1 (red line) and (f}LJo/{F) = 10 (green dashed 
line) with j = 1/2 and Rljq = 1/V0.1 in the gravity-mediated SUSY breaking. 
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FIG. 7: Production rate from the gravity-mediation type of Q ball interacting with the goldstino as a function of 4>oUo/{F) 
for Rljq = 1/V0.2 (red line) and Rluq = 1/V0.1 (green dashed line). The vertical axis is normalized by the saturated rate of 
Eq. flUI- This result can be fitted as (dNi/dt) ~ (6.3+1.6 log(0o^o/(i ? ))-O.O4(log(^>o^o/(f ))) 2 ) X (dN/dt) Eat for i?cj = 1/VTlT. 
The production rates of i = x an d are the same. 

The production rates are shown in Figs. [7] and [8] as a function of 4>o^o/(F) and Rluq, respectively, and they can be 
written as 



log 



6.3 + 1.6 log [ ^ ) -0.04 



w 



tor — 3> 1, /ojq = 



(48) 



10 



We can also get this behavior using Eqs. (j3"5)) and (jUJ) if we naively estimate derivatives in the interaction term as 
d/dr ~ 1/R and ~ oj and use wq > l/R. 
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FIG. 8: Production rate from the gravity-mediation type of Q ball interacting with the goldstino as a function of Rujq 
with <j)oojo/{F) — 0.01. The vertical axis is normalized by the saturated rate of Eq. (|36|l . This result can be fitted as 
(dNi/dt) ~ 1.9 x 10^' s (Rujo — 1.3) x (dN/dt) sa t- The production rates of i = x an d V are the same. 



V. Q BALL DECAY WITH A MASSIVE FERMION 



In this section, we consider the Yukawa interaction with a massive fermion. We are interested in this case because 
squarks interact with quarks and gluinos which have large masses we can not ignore Q. The Lagrangian density is 
written as 



We assign the global U(l) charge for </>, x and r\ such as 1, 1 and 0, respectively. The background field is <j) = 4>( r ) e 
In this section, we take the following background field configuration for simplicity: 



(49) 



— iajQt 



(r) = 4> e(R - r) 



0<r <R, 
R<r. 



(50) 



Inside the Q ball, the Heisenberg equations of motion are 



' ia^dux- qW = 0, 
ia^d^rf - gcf)*x - Mr] = 0, 

ia^d^rj — g<t>x^ — Mrf = 0. 



(51) 



We want to obtain solutions of these equations by the mode and the angular momentum expansion. The symmetry 
of the simultaneous time translations and U(l) rotations allows us to expand the solution as 



X oc e~ lu}t = e~ ik+t , 
X? oc e -*(">- 2w o)* = 



rf oc e 



—i(uj — u}o)t = p — ik v t 



(52) 



where we include x^ and r\ as well as x an d because they mix with each other through the mass term (see 
Eq. (I5ip ). We derive the fermion production rate through the Bogoliubov transformation between creation and 
annihilation operators at t ±oo. Thus, we restrict our attention to < u> < 2wq. 



13 



Now, the solutions inside the Q ball can be written by the spinor solutions in the following form: 

{X = Ae~ lk+t u{k, j, m; r), 
x t = £e ife -*u(/c ! i)77i;r), 
r] = Cer tk ^ t u(k,j, m; r), 
= De~ lkvt u(k,j,m;r), 

where A, B, C, D and k are some constants fixed by solving Eq. (|5ip up to an overall normalization. When we 
substitute these into Eq. (fBTj) . we find 

{(k+ +k)A- g4>D = 0, 
(k v -k)D- g4>A - MC = 0, 
(-k- - k) B - g<t>C = 0, 1 ' 

(k v + k)C~ gcj)B - MD = 0. 

When we eliminate A, B, C and D from these equations, we find that k obeys 

((k + k + )(k - k v ) + g 2 (j)l){{k + fe_)(fc + k v ) + g 2 cf>l) + M 2 (k + k+)(k + fc_) = 0. (55) 

Thus, recalling that fc+, fc_ and k v are fixed in Eq. (|52l) . we obtain four independent solutions inside the Q ball. After 
matching the interior and exterior solutions at r = R, we can write annihilation and creation operators of outgoing 
waves in terms of annihilation and creation operators of incoming waves in the same way as the previous section. 

As mentioned above, we are interested in the case that the scalar field interacts with heavy gluinos. Typically, ujq 
is GeV range in the GMSB, and the mass of gluino is TeV range, soM> u)$. Fortunately, in the case of M > ujq, we 
can calculate the fermion production rate very easily. However, we also present the calculation in the case of M < uq. 



A. Case of M > ujq 



In the case of M > ujq, coefficients T and R of r\ are irrelevant since rj has no degree of freedom outside the Q ball, 
and the boundary condition is rj — > as r — > oo. So, we need to get the coefficients of only \ an d x\ an d the analysis 
is almost equivalent to the previous section. In this case, however, \ an d are related to each other by hermitian 
conjugation. Thus, the Bogoliubov transformation can be written as 



R°Jk 



T°(fc. 



a x (k+J,m) 
(-l)™-4(fe_,j,-m), 



(56) 



/ a x (k + ,j,m) \ 

{{-ly-a^ik-J-m)) ^ v ^ x v .-,„ ^ v .._, jy/ vv ^ ^ 

where R® and T° are coefficients fixed by matching the interior and exterior solutions with the boundary condition 
rj — > as r — > oo, and we write the coefficients in order to maintain consistency with hermitian conjugation. 

We can calculate the \ production rate in the same way as the previous section once we replace ujq with 2wo in 
Eq. (|3"Tj) . The production rate is shown in Fig. |9]as a function of g4>o/uJo, and it can be written as 



dt x 



13 

M 



1.1 x 



dN 
~dt 



for (t? 



'(-) 




\dtj sat 





for 



(-) 

\MJ 



UJ 



— « 1, Rv 



(57) 



> 1, RlUq = TV, 



where we replace ujo with 2wo in the saturated production rate of Eq. (136)) because the energy spectrum of \ i s 
now in the interval (0,2ujq). We can understand the above behavior of the production rate by integrating out the 
heavy particle. Because the effective interaction after integrating out the heavy particle has the form of the Yukawa 
interaction as (g 2 (j)^ 2 /M) (xx/2), the production rate is 12(tt - 1.9)(g 2 ^/Ma;o) 2 x (dN/dt) sat from Eq. J35]) and 
Ruq — n. Thus, we conclude that the production rate calculated in the effective theory is consistent with our 
numerical result of Eq. (|57|) if we replace ujq with 2ojo in the saturated production rate. 

We can explain the above behavior of the production rate in another way. Recall Eq. ([3"5t which is derived in 
the massless case can be interpreted as the decay rate times the charge density times the effective volume, i.e. 
r«£ x (wo<fro) x ^eff- O n the other hand, because the reaction we consider here is two particle scattering process 
4>(p —> qq, the Q ball decay rate should be estimated as (flux) x (cross section) x Wo^q x Ksff- The flux is the number 
density Wo</> 2 times the relative velocity and the cross section can be estimated as g 4 /M 2 from Fig. fTU] Then if we 
assume the relative velocity as O(l) we can get the same parameter dependences with the first line of Eq. (|57|) . 
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FIG. 9: Production rates of x from Q balls as a function of g<f>o/uJo for M/ujq = 2 4 , 2 s and 2 6 with Ruiq = 7r in the Yukawa 
theory with a massive fermion. The vertical axis is normalized by the saturated rate of Eq. (|36[) . 




FIG. 10: Diagram for 



XX- 



B. Case of M < ui 



In the case of M < ljq, r\ can propagate outside the Q ball and its coefficients are also important. Outside the Q 
ball, the fields obey the free equations of motion. The free field expansion of the field \ is t ne same as the previous 
section, but the expansion of the field 77 is complicated by the presence of the mass term. The modes of 77 oc e~ lfe,,t 
and cx e~ lkvt mix with each other by the equations of motion as 



{— ia ■ V77 = — k^rj + Mrf , 
—ia ■ Vr^ = k^ry — Mr}. 

We define the following linear combinations to make 771 and 772 independent from each other: 



_ M fc. ? +P N 

A A 
ky+P M J 



(58) 



(59) 



where p „ 
the solutions are 



hi - M 2 and A = 



Then, with use of Eq. ([22]), 



m = c m inu { 2 >(p,j,m, ;r)e Jfc "* + c mout u {1 >(p, j, m, ; r)e 
m = c V2M u ( - 1 \-p,j,m,;r)e~' tk ^ t + c V2!OUt u^(-p, j, m, ; r)e _lfc "*, 



(60) 



where c Vi (i =1, 2) are arbitrary constants. When \k n \ < Af , p is pure imaginary, and the solutions of 77 and 77^ damp 
outside the Q ball. This is because there is no degree of freedom for 77 at an energy scale below the mass of 77. Thus, 
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we can write the mode expansion as 
r ) = y^ J j dkr, ^c min (k r ,,j,m)e~ %k '> 



A(k n ) 



P=y/k$-M 2 



(terms of c V2 i n ) + (terms of outgoing waves). (61) 



When we quantize the field r], we impose the canonical anticommutation relations for rj and . Then, the coefficients 
c, n and c V2 become operators, and their anticommutation relations are given by {c^. [k rjl j, m) , c ni , (k',j',m'}} = 

8 (k n — fc' ( ) 5jjiS mm >5ii> , where we have used A(k ri ) — <J ^(fc^ + p) 2 + M 2 ^ p/k v and 



M 



%W (p; r)u w {p; r') + (p -> -p) 



dp u wt (p; r)u w (p;r') 



(62) 



The operator c^k^, j,m) is the annihilation operator of the energy E — k rp and its normalization is the same as a x . 

Next, we consider the whole system including the Q ball background. We expand the solution as Eq. ((52j • If k n > 
(k v < 0), we take c m , c V2 {ct, ct) terms in Eq. (|61l) . In the case of < M, rj has no degree of freedom outside the 
Q ball and the analysis is the same as the case of the previous subsection. Thus, the Bogoliubov transformation can 
be written as Eq. (1561) . On the other hand, in the case of | fc^ | > M, there are also incoming and outgoing r\ waves, 
and so the coefficients of rj as well as \ are important. The terms mixing with each other are 



' a xm (k + ,j,m)e- lk + t u { - 1 \-k +1 j,m;r) +(in 
a xm (k^,j,-m)e ik - t (-l) m -u^{-k^,j,m;r) +(in 
c m in(k v ,j, m)e _i Mu( 2 )(p, j, to; r) +(in 



c mm (kr,,j, m)e ^ V 1 ) (_ Pj j ; m; 



+(in 



out, 
out, 
out, 
out, 



« (2) ), 



(63) 



for fc,, > M. After solving the Heisenberg equation of motion and matching the interior and exterior solutions, the 
outgoing waves can be written by the incoming waves. Thus we write the Bogoliubov transformation for the case of 



fc,j > M as 



/ a xout (k + ,m) \ 
(-ir-a^ out (k-,-m) 



V 



L '?72 0Ut 



{k v ,m) ) 



(R xx {k+) T x f x (k+) R mx (k+) R V2X (k + )^ 
T xx f (k-) R x f x i (k-) T mx f (k-) T V2X f (fc_) 
Rxni (k-q) T x \ rn (krj) Rrnvi (k v ) Rq 2 'ni (k-n) 



( a xin (k + ,m) \ 



(-l) m -4in(^-»-^) 
c, 



V 



(64) 



/ 



where we omitted the argument j. The coefficients i?^ and Tj are fixed by matching the interior and exterior solutions 
at 7' = i? and do not depend on m due to rotational invariance. The anticommutation relations imply that the matrix 
of the right hand side is a unitary matrix. 
We can calculate the production rates as 



-^n x (k + ,j,m) = — (0i, 



a lout ( k + > j, ™> XO ut (fc+ , j, m) I Ok 



\T x , x (k + ,j)\ 2 , 

^|T xxt (fc+,j)| 2 + |T I)lX t(fc+,j)| 2 + |T I)2X t(fc+,j)| 2 ) , 



(Af < k^ , 
(|^| <M), 
(k n < —M) , 



(65) 



for the x waves and as 



2tt 



(k v ,j,m) 

0, {k v < M) , 

\T x , m (k n ,j)\ 2 , (Af<fc„<w ), 



(66) 
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for the rj waves. Due to the unitarity of the Bogoliubov transformation, we can see that the production rates satisfy 
the Pauli exclusion principle, < 1 (i = x, Viy Especially, from the second column and the second row of 

Eq. (|64|) . we also have 

|T xx t(fc_,j)| 2 + |T mx t(fc-,j)| 2 + \T V2xf (k_,j)\ 2 = \T xtx (k +lJ )\ 2 + \T x , m (k v , 3 )\ 2 + \T xfr , 2 (k v ,j)\ 2 , (67) 
for M < k n . In other words, using fc + = uj, fc_ = 2uq — u and k n = u — loq, we have 

i=l,2 

where E denotes an energy of each particle and we omitted the argument j and m. From this, we can understand 
the processes of the Q ball decay as the superposition of 

4>{E = uio) + 4>{E = u>o) -> X{E = u) + x(E = 2uj - u), 

0(E = Lj Q )^ X (E = oj)+ m (E = uj o -Lj) 1 (i = 1,2). 1 ' 




FIG. 11: Energy spectrums of the outgoing \ ( re d line), rji (blue dashed line) and r\i (pink dotted line) waves with gfio/^o = 10, 
M/uo — 0.6, Rujq = 7r and j — 1/2. The green dashed line shows dn x /dt(E) — £^ dn Vi / dt(uio — E). 



Fig. [TTI shows an example of the energy spectrum of each field. In this figure, dn x /dt(E) — J^. dn m / dt{u>Q — E) is 
symmetrical; that is, the production rates satisfy Eq. (|68|) . 

Figs. fl~2l and [T3l show the production rates of each field as a function of M/ujq and g</>o/tJo, respectively. We plot 
the x production rate minus the rj production rate instead of the total x production rate; in other words, we plot the 
contribution of only the first line of Eq. (|69[) . From Figs. 1^1 and !12l we can see that the production rate is proportional 
to g 4 (f>Q/M 2 LUQ in the limit of gipo/uo <C 1 and M/luq ^> 1, and it can be written as Eq. (1571) . From Figs. IT^landrTSI 
the production rates in the limit of M /loq <C 1 and g<j)$ /luq 3> 1 can be read as 







2 






X 




























UN 




~ 1.4 x | 







dN 
~~dt 

dN 
~dt 



for^«l, ^»1, i^0=7T. 

OJQ UJQ 



(70) 



17 





0.01 0.1 1 10 



FIG. 13: Production rates of x, an d 772 from Q balls as a function of gcf>o/u)o with Ruo = -n and M/loq = 0.01 in the Yukawa 
theory with a massive fermion. The vertical axis is normalized by the saturated rate of Eq. (|36|) . 



From Figs. IT21 and fT3l the production rates in the limit of M/lo <C 1 and g<f) /uJo <C 1 can be read as 



— ~ 31 

dt x 



UlQ J \ UJQ 



'(-) 




\dt J sat 


cj — >2uj _ 



-^.15[- 



LOQ 



dN 
~~dt 



for M«i, 9jl «i, fla;o=7r . (71) 

LOO LOO 



dt^ * W 



w 



#7. 



We can understand the above behavior of the production rates in the following way. Using the Feynman rules in 
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Fig. 1101 we have the effective interaction, which is a good approximation only in the limit of gipo/tUQ <C 1, as 

M/luo > 1, 
M/luq < 1. 



v 2M 
'g 2 <t>tfM 



XX, 



(72) 



XX, 



Then, substituting this effective coupling (g<p ) x s into Eq. ([38]) gives dN x /dt ~ 12(3 2 ^M/wg) 2 (i?w -1.9) x (dN/dt) sat 
for M/cjq <C 1. This is consistent with our numerical result of dN x /dt in Eq. (|7ip if we replace (dN/dt) sat with 
[(c?A r /c?t) sa t| Wo _ i .2wo]- We also consider the behavior of dN Vi /dt. In the limit of M/ujq <C 1, Eq. (159"]) becomes 



M 
-1 



1 
_M_ 

2w 



(73) 



where we naively take A:^ as the typical energy uiq. From this and Fig. 1141 we have the effective interactions as 



{{gMfxm) ~ [-m/2wo ((aMfxv) + (fctfoffw 1 
{(gMtxm) ~ (-G^o)fx*?) 



(#o)^ ~ Mgfo/uo, 



(74) 



for M/l)q -C 1 and g4>o/ujQ -C 1. Then, substituting these effective couplings (g</>o) into Eq. l[3"8"]). we obtain 
dN m /dt ~ 12( ff( /)oM/a; 2 ) 2 (i?a;o - 1.9) x (dN/dt) sat and dN, n /dt ~ 12(^0 M)) 2 (i?w - 1.9) x (dN/dt) sat . These are 
consistent with our numerical result of dN Vi /dt in Eq. ([71]). Our results indicate that the effective theory gives us 
correct results even for the decay rate of the Q ball. 



x- 



x---- 




FIG. 14: Diagrams for — > x 1 ?- 



VI. APPLICATIONS 



In the GMSB, Q balls can decay only into gravitinos and hadrons if the next LSP mass is larger than luq. We 
apply our results in the previous sections to gravitino and quark production from Q balls in the GMSB and derive 
the branching ratio of the Q ball decay into them. If the gravitino mass is small compared to the typical interaction 
energy, the longitudinal components of the gravitino behaves like the massless goldstino. Since the typical interaction 
energy is uio = O(GeV) 3> m.3/2 in the GMSB, the effective interaction of Eq. (|43|) is a good approximation to calculate 
the gravitino production rate. In this case, we have 

0o^o m s M m 
' c 



(F) 9(F) 

~ (4tt) 2 (F) ■ 

Here we have used Eq. (fTTJ|) and m s M m ~ a(F s )/47r, where a = g 2 /4ir is the fine structure constant, and (F s ) is a 
vacuum expectation value for the F component of a gauge-singlet chiral multiplet in the messenger sector. The SUSY 
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breaking scale (F) needs not be the same as the scale (F s ), i.e. (F s ) < (F). From Eqs. (fTTj) and (]47p . the gravitino 
production rate is calculated as 

Here and hereafter we take m s /ujQ = 10 3 . There is also massive gluino. The Lagrangian density is written as 

£gluino = ^kinetic ~ \/~2g(f)* Xq - -MgXX + h.C. (77) 

where A and q are gluino and quark, respectively, and M g is the gluino mass. Typically, M g — O(TeV) and 
coo = O(GeV), so we can use Eq. (1571) . Using m s ,M m ^> ujq and M g ~ m s in Eq. (ITUl) . we can see that 

(uo/Mg) (V2g4> /uJo) — mlM^/ (Mgul) > 1. Thus, the quark production rate is calculated from Eqs. ((41]) and ((57)) 

as 



d 

^-iV quark ~ 1.1 X 



'(-) 




\ dt /sat 





(78) 



We conclude that the main decay channel is the decay into quarks and is saturated, and the branching ratio of the 
decay into the gravitino is calculated as 



Tfl 

B 3/2 -0.1 (4.8 log— + 7.4 



9 \ f(Fs 



~ 0.08a 

This branching ratio can be rewritten as 



(79) 



TO, \ / TO, Af, 

B 3/2 ~0.1 4.81og— + 7.4 



w o / y \/3g77i3/2-Wp 

to.,M„ 



(80) 



where we use m s M m ~ g 2 (F s ) / (An) 2 and (F) = v / 3m 3/2 Mp (M P = 2.4 x 10 18 GeV: the reduced Planck mass). 

We compare the above branching ratio with the one estimated in Ref. 12]. The quark production rate was 
estimated from the effective coupling g' cS ~ g 2 <f>o/ y / 2irM g ujo for the process squark + squark — > quark + quark via 
gluino exchange. Since we have g' eS (f>Q / (Jq 3> 1, the quark production rate is saturated. The gravitino production 
rate was estimated from the effective coupling g c g ~ u!q/\^2(F) because the elementary process squark — > quark + 
gravitino has the decay rate T = m%/ (\&ir(F) 2 ). Thus, the branching ratio of the decay into the gravitino is estimated 
from Eq. 437]) as 

R , „ gcff^o - r-( (F s )\ 5 m s M m 

B 3/2 - 3?r 0.7-y/a T^r - -p= tt~- (81 

7 wo \\ F ) J y / am 3 / 2 M P 

However, Eq. ((37)) can not be applied to the case of Rujo ~ 1, which is the case in the GMSB, and we should use 
Eqs. ((55]) and fl4"T]) . In addition, we should use (dN / 'dt) S at\u> ^>-2u for the quark production rate because the produced 
quark energy is in the interval (0, 2wo). If we take into account these considerations and use g R, the branching ratio 
is estimated as 

where Eqs. (jSJ) and (|10p are used in the second equality. This result has the same parameter dependences with our 
numerical result of Eq. (|79|) . and the numerical factor is also correct within the order of one. This shows that the 
naive use of the effective coupling is a good approximation even for the decay rate of the Q ball into gravitinos. 
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VII. CONCLUSIONS 



We have calculated the fermion production rates from the step-function type, the gauge-mediation type and the 
gravity- mediation type of Q ball in the Yukawa theory. In the limit of g4>o/ 'u>o < 1, we have found that the decay rates 
can be obtained from the step-function type of Q ball with the correction factors coming from the relations between 
the total charge Q and the radius R ~ 1/u)q. On the other hand, the decay rates in the limit of g<fio/oJo ^> 1 are 
saturated and proportional to the square of the effective Q ball radius R! , where R' is determined by g4>(R') / uj ~ 1. 

We have also calculated the goldstino production rates from the gauge-mediation type and the gravity-mediation 
type of Q ball using the low energy interaction with the supercurrent. Our results can be explained by the pro- 
duction rate through the Yukawa interaction where the Yukawa coupling g is replaced by the effective coupling 
g c s ~ ujq/\/2(F). This effective coupling comes from the fact that the elementary process squark — > quark + 
gravitino has the decay rate T = m^/ (16ir(F) 2 ) or we can naively estimate derivatives in the interaction term as 
d/dr ~ 1/R and d/dt ~ w . 

We have also calculated the Q ball decay rates in the Yukawa theory with a massive fermion. Our results are 
consistent with the effective theory once we make the replacement u>q —> 2u>o in the saturated rate of the massless 
fermion, since the produced fermion energy is in the interval (0, 2u>o). Especially, when ujq < M, we can integrate 
out the heavy particle and use the effective coupling {g4>o) cS — g 2 (j) 2 ] /2M in the Yukawa theory with the replacement 
ujq — > 2wo in the saturated rate. 

In the GMSB model, the branching ratio of the decay into the gravitino has been calculated as 



for m s /ujq = 10 3 . This branching ratio is much less than the one estimated in Ref. \i2\ - The main reason is that 
Ref. 12] used the production rate for the limit of Rw — > oo which is not valid in the GMSB. Another reason is that 
we should take into account the Q ball configuration which is different from a step function. Therefore, the gravitino 
dark matter from the Q balls in the GMSB should be reconsider using the correct decay rates obtained in the present 
paper, which will be presented elsewhere (l9j . 



We thank Shinta Kasuya for useful discussions. This work is supported by Grant-in- Aid for Scientific research from 
the Ministry of Education, Science, Sports, and Culture (MEXT), Japan, No. 14102004 (M.K.), No. 21111006 (M.K.) 
and also by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan. 



[1] I. Affleck and M. Dine, Nucl. Phys. B 249, 361 (1985). 

[2] M. Dine, L. Randall and S. D. Thomas, Nucl. Phys. B 458, 291 (1996). 

[3] A. Kusenko and M. E. Shaposhnikov, Phys. Lett. B 418, 46 (1998). 

[4] K. Enqvist and J. McDonald, Phys. Lett. B 425, 309 (1998); Nucl. Phys. B 538, 321 (1999). 

[5] S. Kasuya and M. Kawasaki, Phys. Rev. D 61, 041301 (R) (2000). 

[6] S. Kasuya and M. Kawasaki, Phys. Rev. D 62, 023512 (2000). 

[7] S. Kasuya and M. Kawasaki, Phys. Rev. D 64, 123515 (2001). 

[8] A. Kusenko, L. Loveridge and M. Shaposhnikov, Phys. Rev. D 72, 025015 (2005) 

[9] A. Kusenko, L. C. Loveridge and M. Shaposhnikov, JCAP 0508, 011 (2005) 

[10] I. M. Shoemaker and A. Kusenko, Phys. Rev. D 80, 075021 (2009). 

[11] F. Doddato and J. McDonald, JCAP 1106, 008 (2011). 

[12] S. Kasuya and M. Kawasaki, Phys. Rev. D 84, 123528 (2011). 

[13] A. G. Cohen, S. R. Coleman, H. Georgi and A. Manohar, Nucl. Phys. B 272, 301 (1986). 

[14] J. Hisano, M. M. Nojiri and N. Okada, Phys. Rev. D 64, 023511 (2001). 

[15] S. Coleman, Nucl. Phys. B262 (1985) 263. 




(83) 



Acknowledgment 



21 



[16] A. Kusenko, Phys. Lett. B405 (1997) 108. 

[17] A. de Gouvea, T. Moroi and H. Murayama, Phys. Rev. D 56, 1281 (1997). 
[18] T. Multamaki and I. Vilja, Nucl. Phys. B 574, 130 (2000). 
[19] S. Kasuya, M. Kawasaki and M. Yamada, in preparation. 



